# Exercicio 8 - Regresses Multiplas
# Camila Galheigo Coelho 

setwd("C:/Users/Camila/Documents/PESQUISA/1_DOUTORADO/Disciplinas/R")
getwd()

# Galileu estava certo? #
init.h = c(600, 700, 800, 950, 1100, 1300, 1500)
h.d = c(253, 337, 395, 451, 495, 534, 573)

modelo.linear<-lm(h.d~init.h) 
plot(h.d~init.h)
abline(modelo.linear)

modelo.2o.grau<-lm(h.d~init.h+I(init.h^2)) #2o grau
curve(coef(modelo.2o.grau)[1]+coef(modelo.2o.grau)[2]*x+coef(modelo.2o.grau)[3]*x^2, add=T, lty=2)

modelo.3o.grau<-lm(h.d~init.h+I(init.h^2)+I(init.h^3)) #3o grau
curve(oef(modelo.3o.grau)[1]+oef(modelo.3o.grau)[2]*x+oef(modelo.3o.grau)[3]*x^2+oef(modelo.3o.grau)[4]*x^3, add=T, lty=3)

summary(modelo.linear)
summary(modelo.2o.grau)
summary(modelo.3o.grau)

anova(modelo.linear,modelo.2o.grau,modelo.3o.grau) 

#a menor soma dos quadrados dos resduos  a do modelo de 3o grau e portanto este  o que melhor.


# Recem nascidos #

RN<-read.table("babies.txt",head=TRUE, sep="\t")
head(RN)

save.image()